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ABSTRACT 


A five degree-of-freedom linear model for the XR-3 
surface effect ship is developed for constant speed operation. 
The weight removal transient response and the roll behavior 
of the craft in turn maneuvers are investigated by both 
linear and six degree-of-freedom nonlinear models, and the 
linear model simulation results are compared with the non- 
linear model simulations. Some of the parameters of the 


craft are investigated for roll motion. 
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I. INTRODUCTION 


A. BACKGROUND 

The Surface Effect hybrid type high performance ship has 
been of great interest recently, and in 1970 the XR-3 test 
craft was delivered to the Naval Postgraduate School. As 
seen in [7] Leo and Boncal converted the digital simulation 
program for the 100-B test craft to represent the XR-3, thus 
providing a nonlinear model. Іп [1], Gerba and Thaler have 
developed a heave-only linear model and further studies were 
made by Barnes in (2] by developing a two degree-of-freedom 


linear pitch and heave model representing the XR-3 test craft. 


В. OBJECTIVE 

In order to investigate the roll behavior of the XR-3 
test craft during turn maneuvers, a five degree-of-freedom 
linear model is developed in this thesis. Some modifications 
are made for pitch and heave modes, therefore test of the 
vertical plane by means of weight removal is reinvestigated 
and compared with the nonlinear simulation results. 

In [3], the roll behavior of the XR-3 and the sidewall 
sensitivity are studied using a six degree-of-freedom non- 
linear model. Similar studies are made in this thesis after 
developing a five degree-of-freedom linear model and comparing 
with the nonlinear model results. Finally, some modifications 


are suggested for nonlinear modeling. 
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II. NONLINEAR XR-3 MODELING FOR FIVE DEGREE-OF-FREEDOM 


A. ASSUMPTIONS 

The following assumptions and simplifications are used 
during development of a constant speed XR-3 test craft model: 

1. Aerodynamic drag on the forward seal, and superstruc- 
ture skin friction forces are neglected. 

2. Momentum (ram) drag is disregarded for all modes. 

In other words, it is accepted that the location of the lift 
fans with respect to the craft C.G. could not affect signi- 
ficantly the roll and pitch stability. Also, air stagnation 
force in the z-direction is assumed to be negligible. 

3. During dynamic response of the craft, only air 
leakage with constant leakage area at the stern seal is 
assumed. No other air escape is considered such as possible 
leakage from sidewalls or bow seal. 

4. Plenum pressure and air mass thermodynamics are 
represented by ideal diabatic processes with the air constant 
specific heat ratio y = 1.4. 

5. Effective plenum roof area is developed as a function 
of craft weight. 

6. Planing forces are concentrated at the centroid of 
the keel longitudinal cross section through the fluid medium. 

7. Pitch and roll damping moments due to added mass are 
included in the craft dynamics, but freue contribution to 


heave motion is neglected. 
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8. The thrust forces and all hydrodynamic resistive 
forces (drag) in the z-direction cancel each other, since 
at constant speed operation, these forces are in equilibrium. 
However they might contribute to the moment equations such 
as pitch moment, but considering their effective moment arms 
with respect to the C.G. are almost equal to the thrust moment 
arm, so that their contribution in the moment equations is 
also neglected. 

9. For seal modeling, vertical resistive forces due to 
hydrodynamic response are taken into account including 
corresponding effective wetted seal area. 

10. Simplified craft geometry is used for plenum, 


buoyancy and seal modeling. 


B. COORDINATE SYSTEM AND SIGN CONVENTIONS 

1. The cartesian coordinate system is used with its 
Origin located on the water line corresponding to the craft 
center of gravity. 

2. The x-axis lies on the water line with the positive 
longitudinal displacement being measured forward. 

3. The y-axis lies in the vertical cross section plane 
with the lateral dislacement being measured positive to 
starboard. 

4. The z-axis lies in the vertical cross section plane 
with the positive vertical displacement being measured 


downward. 
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5. Z represents heave C.G. displacement measured from 
the origin (at water line) to the center of gravity. 

6. Z 15 one of the craft dimensional values measured 
from the keel to the C.G. 

7. Pitch is the rotational motion about the y-axis, 
angular pitch displacement is represented by 0, positive 
pitch angle is measured upward (boat noses up). 

8. Roll is the rotational motion about the x-axis, 
$ angle represents angular roll displacement, positive roll 
angle is measured to starboard (boat heels to starboard). 

9. Similarly the yaw is the rotational motion about the 
Z-axis, positive yaw is measured in the clockwise direction. 

10. Positive roll, pitch and yaw moments are to be 
defined as the moments in the positive roll, pitch and 
yaw directions correspondingly. 

ll. As a reference zero pitch and roll angle provide the 
x-y plane parallel to the water surface. 


The craft coordinate system is shown in figure 1. 


С. SIMPLIFIED CRAFT MODELING AND GEOMETRY 


1. General Simplified Geometry and Dimensions of the 
XR-3 Test Craft 


In figure 1 the sideview simplified geometry of the 
XR-3 test craft is shown. Also in figure 3 the top view is 
illustrated. The variables and dimensions in these figures 


are defined as follows; 
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Figure 2. Bow seal geometry 


Aft section Forward section 





TOP VIEW 


Figure 3. Simplified geometry (top view) 
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& = Draft of the craft at center of gravity 


C = Plenum pressure center (the centroid of the 
effective plenum area) 


Rorke = The longitudinal distance of the sections’ 
centroids from the C.G. (forward and aft 
section centroids respectively) 

Lara = Bow and stern seals front face distances 
at water level 

Lo = Centroid distance of the planing force 
effective area 

2,12, = Forward and aft sections sidewall submerged 


region length. 

All craft dimensional values are taken from the XR-3 
six degree-of-freedom model scale drawings as indicated in 
[2] and [3]. 

2. Sidewalls 

As reported in [2], the sidewall modeling is developed 
from the SIDEWALL subroutine of the 6 DOF computer model. 
Like craft dimensions, all sidewall dimensions and their 
related geometric configurations are based on scale model 
drawings. For simplicity of calculation, the sidewalls are 
divided into two sections called forward and aft. During 
dynamic response of the sidewalls, these sections are con- 
sidered individually with their own centroids, so related 
calculations can be done easily. Especially, the sidewall 
modeling is developed in order to determine buoyancy forces 
and moments. By means of this modeling, the average sidewall 


width, average draft at section centroid and effective 
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plenum area are also determined and used extensively in the 
rest of the calculations. 

The sidewall vertical cross sections at the section 
centroid forward and aft are sketched in figures 4(a) and 4(b) 
respectively. From these figures, the average forward side- 


wall width at the buoyancy center is found to be 


x 


ES f 
b. EN 5. + b. (II-1) 
where 
Я = Draft at forward buoyancy center 
tan B, =  Tangent of the sidewall deadrise angle 
at the forward buoyancy center 
De = Sidewall flat surface width at the kell 


(forward section) 


Therefore submerged sidewall volume of the forward 


section is given as follows 
Ve er (II-2) 
ὅς can be formulated as 


5, =ош = апо 


then the small angle approximation tan 0 = 0 gives 
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Forward section Aft section 


Figure 4(a) 





Figure 5. 


Figure 4(b) 


Plenum modeling 


19 


Vw... 





9 (II-3) 
As defined before, La draft at C.G. is 


ГА = 2 +2 (II-4) 


Similar expressions can be derived for the aft 


section, 


ор πι ο 


where 


5. = sidewall flat surface width at the keel 
(aft section) 
= E 
У, δ. 8» 8. (II-6) 
A = Ха + Xe 8 


Note that Ler 5. are functions of draft and pitch 


angle. 
3. Plenum Modeling 

In [2] plenum modeling is presented including plenum 
pressure variation due to changing of draft. In the develop- 
ment of a pitch, heave and roll linear model, the simplified 
plenum geometry and dynamics which were derived by Barnes 
are considered and used in our plenum dynamics representation. 

The variables and dimensions for plenum modeling, 


as seen in figure 5, are found as 
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= Е -ᾱ E 
NS = Ro Ro (=) (II-8) 
т b. КЕ (11-9) 
59 ü MB = 5 
where 

νο = Length of plenum at the water line 

А, = Effective plenum area 

b = Lateral distance from one sidewall to 


another in the plenum region. 


ο Lor 5а! о and h are all constants related to 
craft dimensions. Only the effective plenum area A, varies 
when the craft draft varies, that means variation of the plenum 
pressure centroid Co is forward or backward. 
4. Bow and Aft Seals 

The seals are an involved additional dynamic system 
with their own degree of freedom relative to the craft, and 
they carry little load compared with the other components 
of the craft. However in the simplified model, the seal 
motions relative to the craft are neglected. 

The assumed seal geometry is shown in figure 2, 
developed in [2]. Using an appropriate geometric configura- 
tion, the seal wetted area is derived as a function of craft 


draft and pitch angle. 
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From [2], the width of the wetted surface for the 


bow seal is; 


£. ~ £a tan Ө 
_ ^d 7 ^3 : 
Εως sin 31° (11-10) 


where 


sin of the angle between water surface 
and the bow seal front face. 


sin 31° 


The numerator of this ratio represents draft at the 
bow seal front face. 


A similar expression for the stern seal is; 


Х. +4, tan Ө 


ο а. 
а sin 32° 


Multiplying ХЕ апа X, by bu bow and stern seals 


wetted areas are obtained as follows; 


5. Rudder and Propulsion Hydrodynamics 


The rudder hydrodynamics is calculated from [4], 
inserting rudder geometry dimensional values into the fin 
equations. Actually we are interested in the lateral force 
(sway) on the rudders for a given rudder order deflection. 


This dynamic equation is given in [4] as follows; 
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Y = 
r 


where 


2 
p u A (1 + 


d 
ё 


h 


um Уһ 
) C ló, = (1 B) zu (11-11 (8) ) 
I 12 


Rudder platform area (for one rudder only) 


draft at the stern = z + 2. _ X, 8 


Vertical depth to bottom tip of the rudder 


Rudder lift coefficient 


27 R 
a 


R + 3 (for fully wetted flow regime) 
a 
The rudder aspect ratio 


Rudder deflection (the sign convention of the 
rudder is to be positive for port turn) 


MISC Tm р 


r S 


longitudinal coordinate of rudder geometric 
center (always negative value) 


lateral speed (ft/sec) 


yaw rate (rad/sec) 


In order to simplify the rudder dynamic equation, 


the following coefficients are defined as 


Бя 
ат 
г 
= ιν 
DM E r η. Υ 
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Then equation (II-ll(a)) becomes 


Mi 
Y. = C. (ó, = Се =) (II-11 (b) ) 
Let us define 2. to be the vertical coordinate of 
the rudder from the C.G. to its geometric center, the rudder 


roll and yaw moments can be calculated from 
K = = Y 7 (ІІ-11(с)) 
М = Y x (11-11 (а)) 


During turn maneuvers, the propulsion thrust vector 
is also deflected according to the given rudder order. There- 
fore the lateral component of the thrust contributes yaw and 
roll moments as well as a sway mode, but this diturbance must 
be considered as a step input to the linear modes with constant 


magnitude, therefore the total propulsion sway force is 


Yor = 2 X, sin ὃ, (II-12) 


where X, is the propulsion thrust corresponding to craft speed. 


Similarly the roll and yaw moments generated by the 


propulsion system are 


K = У 2 (II-13(a)) 
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N = Y x (II-13(b)) 


where < and Zo are the longitudinal and vertical locations 


of the propeller respectively. 


D. DYNAMIC FORCES FOR FIVE DEGREE-OF-FREEDOM MODES 
l. Seal Forces 
Considering the seal configuration and dynamics, the 


pressure difference across the wetted portion of the seal is 


Б ол p U ta 
where 
p = Plenum gauge pressure (1b/£t^) 
P, = Atmospheric pressure 


Since the seal is in equilibrium, the vertical hydro- 


dynamic force on the seal is 


sf b £f 
$. - 2.08 
Nem, d а” 73 Е 
EE ba Csin που (II-14) 


For the stern seal hydrodynamics, a similar expression 


can be found as 


B 84 +240 
а а aa 


sa (II-15) 
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where 


P. = Additional pressure supplied into stern air 
bags in order to reduce leakage. 


2. Buoyant Forces 


The forward section buoyant force is determined from 


equation (II-2), it is 


2 РУ 


bf 1 


here p is the density of the seawater and g is the gravitational 
constant. 


Finally, the forward and aft section forces are 


А 
TA ε o. 3 
δε Ξ 2 р 9 (5 tan 6, + be) A. Ly (II-16) 
La 
а "© 2 95 вап В. + bafa 52 (II-17) 


3. Planing Forces 
In addition to the pitch angle 0, the sidewall flat 
area at the keel determines the planing effect of the craft. 
When the craft pitches up, the frontal area (through the fluid 
medium) causes a lift-drag force in the z-direction. 


The planing force, from Reidel's work [3], is; 


- Д x 
207, = > T p Au sin Ө 
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where 


flat surface area at the keel 


p 
|| 


speed of the craft 


с 
|| 


Using the small angle approximation and including the 
loss coefficient Col suggested by Barnes in [2], the planing 


force is redetermined as 


2 = = CÇ 


pl pl το u б (II-18) 


A similar expression is derived for sidewall dynamics 

in [4], but in our modelling equation (II-18) is used. 
4. Pitch Damping Forces 

The sidewalls of the surface effect ship influence 
moderately heave, pitch and roll modes, so in [4] the hydro- 
dynamic and hydrostatic forces and moments acting on the side- 
walls of uniform triangular cross section are derived with 
particular application to vertical plane motion and lateral 
motion. 

Considering constant speed operation of the craft, 
neglecting some terms in the equations provided from [4], 


the vertical force acting on only one sidewall is found as; 


(II-19) 
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where 


А. = vertical added mass calculated at the stern 
2 
- (7/8) p B, 
B. = value of sidewall waterline at the sidewall 
stern 
X. = stern distance from C.G. (negative value) 
а = pitch rate (rad/sec). 


For the pitch damping moment, the same equation is 
used by Barnes in [2]. 
5. Roll Damping Forces 
The XR-3 test craft is dynamically similar to 100-B 
surface effect ship. Therefore we can use dynamic equations 
calculated in [4] for the XR-3 test craft. The heave-roll 


damping force for one sidewall can be written as; 
7 =. - A... u b_ p (II-20) 
In the above equation, b. is half of the craft width 
and it must be used (+) for starboard sidewall, (-) for port 
sidewall. 
Similarly the lateral damping force due to sidewall 


hydrodynamics is 


Y = -A u(v + x r _ z P) (II-21) 
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here A), is the lateral added mass calculated at the stern, 
p is the roll rate. A), is calculated from [4], that is 
A = 1/5 ра В 
22 5 
6. Plenum Pressure Lift Force 
From [2] and [4], the plenum pressure lift force is 
determined assuming uniform pressure distribution, This is 


formulated as 


Z, = A Pb, 
from equations (II-8) and (II-9) follow 


E & Ra 
2 = -= b. P (Lo = -2-) (II-22) 
E. DYNAMIC PLENUM AIR MASS, PRESSURE EQUATIONS 

Assuming constant aft seal leakage area, the leakage flow 


rate (Q ) is given by Gerba and Thaler for the XR-3 heave 


out 


only model in [1], then we have 


2 P 
2 b, 1/2 _ 


where 


Nozzle flow coefficinet 


а 
N 


> 
u 


Leakage area 
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DET = Atmospheric air density. 


The linearized air flow (O μ) into the plenum is pre- 
sented in [5] according to the fan properties and calculated 


in [1] as follows, 
(II-24) 


where n is the number of fans supplying the air to the plenun, 
B. is the steady state air flow rate corresponding to steady 
gauge pressure, k. is (1/з1оре) of the РЬ VS. Q;, curve at 
the equilibrium point. 


From equations (II-23) and (II-24), the plenum air flow 


rate can be written as 


. dm, 
= Pa Qin T Cont? 115) 


The equation (II-25) reveals thatif the fan flow rate 
(Q; y) equals to the leakage flow rate er then the plenum 
air mass (M, ) is constant. 

The absolute plenum pressure is calculated from 


Mp 


- 
P = P (L——) (II-26) 
b a Vo Pa 


In the above equation, УЬ is defined as plenum volume 


and its simplified form is calculated from [1], 
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M = ou E A 5а! (11-27) 


where Va is to be the empty plenum volume. 


F. NET HEAVE FORCE 
The resultant heave force is the summation of the all 


external dynamic forces, it is 


where 


Z sw also includes a pair of heave buoyancy forces due to 
roll, but they are in opposite directions. Therefore these 
stabilizing forces are only taken into account in the roll 


motion. 


G. PITCH AND ROLL MOMENTS 
l. Seal Moments 
a. Pitch Mode 
The seal force is modeled by Barnes in [2] acting 
on the center of the seal wetted length. The effective pitch 


moment arms are;j 


Log = 2 _ хе /0 (for forward seal) (ТТ-28) 


> 
II 


E T д + х_/2 (for aft seal) (II-29) 


S 








From this follows the forward and aft seal moments 


Mor = = Los er (II-30) 


М a = -= zr S (1I=31) 
where Dos and 255 have already been defined. 
b. Roll Mode 

The flexible seals can contribute little to roll 
stability. Actually the seal configurations and hydrodynamics 
are relatively more complex. In [4], the seal roll moments 
are calculated by combining the roll contribution of all seal 
elements. 

Considering the seal response, the submerged 
region of the seal can be calculated and its roll moment con- 


tribution can be formulated in the following equations; 


3 
= b. и вах 


A 
II 


tan ф (11230) 


sf f 


3 
- b, /12 ра x, tan ф (II-31) 


A 
li 


Sa 


For this modeling, the roll moment arm is considered 
as (2/3 б), that is the centroid of the triangle cross section 
of the submerged portion. 

2. Buoyant Moments 
a. Pitch Mode 
The pitch buoyancy moments are calculated from 


equation (II-16) and (II-17) assuming constant moment arms 
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T and 16: The pitch buoyancy moments are 


M E = - Ls d (II-32) 


ll 
са 
> 


Mpa ba 5 (11-33) 


b. Roll Mode 
As discussed in the heave mode, the effect of the 
sidewall buoyancy due to the roll is to be a stabilizing moment. 
From figures 4(a) and 4(b) the buoyancy volume change can 
be approximated as 
be tb. A Q 


f f 


au EE 2 ш 2 tan 8? can В, 


£ +65.) 2$ (11-34) 


Finally, roll moment is 


K = — 2 p g AV 


bf (II=35) 


£ 


Similarly, the aft section roll moment is 


К м 204 AV, 


3. Planing Pitch Moment 
The planing effect due to the sidewall flat area at 


the keel contributes only pitch motion as follows; 
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4. Pitch Damping Moment 


This is calculated from equation (II-19) as 


In the above equation, the factor "2" appears because 
there are two sidewall sections. 
5. Roll Damping Moment 
Combining equation (II-20) and (II-21), the roll 


damping component is 


K = 27 D. - У z +K (II-37) 
s wave - 


where Ка represents the additional roll damping term 
because of vertical wave generation in the roll. 


From [3] and [4], Кс 15 


K ое в в (II-38) 
wave О C 


Here, B, is the sidewall vertical added mass given as 


Вс = π/δ ρ 48 4 2 ο š 


da 
6. Plenum Pitch Moment 
From equation (II-22) and defining Хор to be the 
plenum pressure centroid, the plenum pitch moment is given 


as 
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M = — Z X = A, (P 


p р==ср b 7 Pa? А 


b p 


where X. is calculated as a function of craft draft in [2] 


P 
as follows; 


Хор = X. + (2 & )/2h 


X, л Х.д. 


7. Plenum Roll Moment 
The destabilizing static plenum roll moment is 
generated when the craft rolls to one side. Since the effec- 
tive plenum pressure centroid shifts laterally to the other 
side tending to increase rolling of the craft. The plenum 


roll moment is 
Ko = A pPp - Pi) (z, - 54209 
8. Rudder and Propulsion Roll Moments 


These roll moments have already been calculated 


from equation (II-11(c)) and (II-13(a)). 


Н. SWAY FORCES 
1. Cross Flow Drag (Sidewall) 
During the turn maneuver, the cross flow drag is 
generated acting on the outboard sidewall. The direction 


of this force is negative while turning to the left. 
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Neglecting the lateral drag component inside of the sidewall, 


the net cross flow drag can be calculated from 


ПЕР І. 2 _ 
D ~” 5 р sign (V_) Са V. & (£a Е b $) (II-39) 
where 
МЕ ЫК rus (2. = &4/2)p 
Са = cross flow drag coefficient 


Considering the outboard surface of the submerged 
sidewall sections (dead rise angle effect),the lift force 
in the z-direction can be produced due to the lateral drag 
deadrise projection. In [4], the importance of the sidewall 
deadrise angle is mentioned, however at the hydrodynamic 
section in this reference, the lift component of the drag 
force has not been derived. 

In this modeling, the lift force (24) and lateral 


drag component is calculated assuming both have the same drag 


coefficients (Са = 1.28). They are 
Ya = D sin B (II-40) 
24 = Dau; cos В | (II-41) 
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From equation (II-40) and (II-21), the total hydro- 


dynamic sway force of the sidewall is 


Y ow = (D w + GA sln B+ ‘sw, (II-42) 


2, Centrifugal Force 


The centrifugal force acting on the C.G. of the craft 
is always in the outboard direction during the turn maneuver. 


Ус is the centrifugal force, given as 


u? 
лт ле 


С 
R, the turning radius is equal to (u/r), combining this 


with the above equation gives 
Y = -mur (II-43) 


a negative sign is used due to the coordinate system sign 
convention. Since the right turn produces (+) yaw rate, then 
Yo is calculated as (-), so outboard sway force is generated 
for this configuration. 
3. Rudder and Propulsion Sway Forces 

They are determined from equation (II-11(b)) and 
(II-12). 

The sway forces due to the seals, air plenum and 


aerodynamic drag are neglected in this modeling. 
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I. YAW MOMENTS 


1. Sidewall Yaw Moment 


The sidewall yaw dynamic equations are given in [4]. 


After some simplification, the net widewall yaw moment is 


calculated from 


Now = N, + N, + (Ko = Χ Ὁ) bo 
where 

N, = 2Zu(v+x, 7-2 PJA,nA(l-A) 

OO υπ р) 
d ο 5 
E 3 2 2 
N. = О Са Ха (2 Si $9 r V + 5, L V. ) 
Q + X. 

а = мл СС The ratio of C.G. location of 
sidewall bow to the craft 
length 

E DEAN DOE» 

1 3 
| 2À - 1 
πω; 
X = - 1/2 С.А ue 
pw P “£ “q 
p 
= port wetted area surge resistive force 
X = — 1⁄2 p C- À 42 
sw Е We 


= starbord wetted area surge resistive 
force 
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(II-44 (a) ) 


(II-44 (b)) 


(ІІ-44(с)) 


(ІІ-44(а)) 


(ІІ-44(е)) 


(II-44(f)) 


(II-44(g)) 


(II-44(h)) 





a, = [2 (4, £ b $) + b 1k (II-45) 
b. +b 
2 f d 
bs E 2 


Average sidewall flat surface width at the keel. 


In equation (II-45) * term related to the starboard 
or port wetted surface calculations. 

2. Rudder and Propulsion Yaw Moments 

The yaw moments due to the rudder and propulsion sys- 
tem are determined from equation (II-ll(d)) and (II-13(b)) 
respectively. 

The seals, air plenum and aerodynamic yaw moments are 
neglected. Actually the aerodynamic yaw moment should be taken 
into account. Since during the turn maneuver, the yaw drift 
angle is produced, thus a positive aerodynamic yaw moment 
contribution to the yaw motion has been generated. From 6 
DOF nonlinear model simulations, this contribution is determined 


at 5.2% of the total yaw moment. 


J. TOTAL PITCH MOMENT 
The summation of the all dynamic pitch moments gives the 


net pitch moment which can be written as 


M = M_ £ + e t M) £ + Mba + sib 4 М, + E (II-46) 
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K. TOTAL ROLL MOMENT 

Before calculating net roll moment, the roll moment 
equation due to the cross flow drag term must be derived. 
Including also the added mass effect to the equation (II-40) 
and (II-41), the result is 


Ka = ο + Ysw? Po cos ß - Zo sin B) (II-47) 


Combining all roll dynamic equations, the net roll 


moment is calculated as 


K = K + К + К + K + К + K т K, + Кр К 


sf sa bf ba sw pr а 


L. TOTAL SWAY FORCE 
From equation (II-ll(a)), (II-12), (тт-42) and (II-43), 


the total sway force is obtained as 
ү = ү +ү + Y + Y (II-48) 
M. TOTAL YAW MOMENT 
The summation of the sidewall, rudder and propulsion 


yaw moments gives the net turning moment about the z-axis 


N = Now + N. + Nor (II-49) 
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N. FIVE DEGREE-OF-FREEDOM EQUATIONS OF MOTION 

Considering equations of motion for XR-3 test craft, 
the yaw (1) is defined as the rotational motion of the craft 
coordinate system relative to the fixed frame and the positive 
U angle is measured in the clockwise direction. This sign 
convention is in agreement with the rudder order sign defini- 
tion. Since positive rudder angle corresponds to turning 
to the left, producing negative yaw moment and finally at 
the equilibrium condition, the negative yaw moment is 
generated. 

Assuming small values of pitch (0) and roll angles and 
relating heave motion to the free surface, the relationship 
between the fixed and craft coordinate system can be written. 
Eliminating the terms related to the fixed coordinate system, 
the simplified 5 DOF equations of motion are obtained as 


follows; 


mv - Y (II-50) 

m м = 7 (11-51) 

Ге. р EI. r = K (11-52) 
уу q = M (ТТ-53) 

I. r - I,, p = N (II-54) 
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The right hand sides of the equations of motion give the 
forces and moments acting on the craft in terms of the craft 
in terms of the craft coordinate system variables. In 
equation (II-57) 75 represents heave acceleration and all 
other variables in these equations have already been defined. 

Im order to develop a 5 DOF linear model, only the surge 
equation of motion (x-direction) is neglected. Since from 
6 DOF nonlinear simulation results, very small velocity 
reduction is observed during the turn maneuver of the XR-3 


test craft (over the speed range 20 through 30 knots). 
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ITI. LINEARIZATION OF THE EQUATIONS 


A. TAYLOR SERIES EXPANSION AND TOTAL DIFFERENTIAL 

Generally, the linearization process can be performed by 
means of a Taylor series expansion using the linear part of 
the series (omitting high order terms). The same process can 
be easily done applying a total differential to the net force 
and moment equations, since the effect of small disturbances 
to the system is to be determined. 

The total differential of a given dependent variable A 
with respect to the independent variables x, y and z is given 
as 
9A 


dy + 35 dz (ттт=1) 


dA = ЗА ax + 
9 x ду 


B.  LINEARIZATION OF THE EQUATIONS OF MOTION AND AIR FLOW 
DYNAMICS 


The total differentials of the equations of motion and 


air flow dynamics have been derived as 


dw = dZ/m (III-2) 
dq = AMI y (aras) 
. p Е 
dm, = -ipa P Kp + C, л, (——®—)1/?уару Ir 
2P, (0) 
dv = аут | (Tira) 
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where Р (0) represents the steady state value of the plenum 
gauge pressure according to the final equilibrium condition. 

The equations (II-52) and (II-54) both involved second 
order yaw and roll acceleration terms. “herefore applying 
the total differential and solving these two equations 


for dp and dr we have 


Ω, 
e 
li 


dK/I + dN/I, (III-6) 


Qu 
rO 
| 


ак/ту + ам/т (III-7) 


E, I, and I. are defined as 


I 
XX ZZ XZ 


li> 


(III-8 (a) ) 
XZ 


I 
_XX 22 Х2. (III-8 (b)) 
xx 


H 
i> 


O т 2 
т, 5 HA (III-8(c)) 


22 


Equation (III-6) and equation (III-7) contain both net 
incremental roll (dK) and yaw (dN) moments. This is simply 
a cross product inertial term (1.2) effect іп the equations 
of motion. 

In equation (II-59), ар, is not a state variable, it сап 


be considered as an output variable, therefore from equation 


(II-26) and (II-27), ар, must be 
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C. FIVE DEGREE-OF-FREEDOM LINEAR MODEL SENSITIVITY COEFFICIENTS 

The sensitivity coefficients related to the right hand side 
of the equations of motion are defined as the partial deriva- 
tives of the net forces and moments divided by craft mass or 
inertial terms. 

In [2], Barnes has developed a 2 DOF (heave and pitch) 
linear XR-3 model. Therefore most of the neave and pitch 
sensitivity coefficients are directly used in order to develop 
a 5 DOF linear model. However some modification is made for 
plenum sensitivity coefficients and pitch mode. Also some 
of the steady state operating point values have been reinves- 
tigated using the 6 DOF nonlinear model. 

In all the linearized equations, the final steady state 
values of the variables have been used. 

1. Pitch and Heave Sensitivity Coefficients 

In [2] equation (III-40), the sensitivity coefficient 


due to z has been redefined as 


ο. - b, Lo 
DHPZ = - AL (0) s= + PL (0) h 





This is claculated from (II-22) where the second term 
does not appear in [2]. 
The modified form of the plenum pitch sensitivity 


coefficient due to state variable z is to be 
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dg — by Х9 
DPPBZ = - 2,10) 2R ^ а Р, (0) Б 





The derivation of these modified equations is based 
upon the total differential of the plenum lift force, it 


15 


92 92 
az = —Рар + а 
p ӘБ, b ЗЕ D 


From equations (II-9), (II-26) and (II-27), the total 


differentials ау, апа ар, have been calculated as 


Do. Йй 
22 73 Е м 9 Z 
dV, = (A, (0) & (0) Б ) dz (III-9) 





= Ж 1 w 9 
dP, - ή Pp (0) I (gy d: үш s το. Б )142. (ІІІ-10) 





ар, апа av, equations are also different from the equations 
derived in [2]. 

The bow and aft seal pitch moments sensitivity 
coefficients (DPSFTH and DPSATH) due to 6 and plenum air 
mass equation (AMB) are rederived and their modified forms 
are used for the pitch mode in the 5 DOF linear model computer 
program. 

Considering the 5 DOF linear model state space matrix 
representation, the ajj S represent "A" matrix elements. 


Therefore the ajj values relating to the sensitivity coeffi- 


cients for heave motion are determined as 
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821 = 0225 
Άρη = DZTH 
а = DZDP 
а,5 = DZMB 
where DZZ, DZTH, etc., have been derived in [2]. However 


their modified forms can be found in the 5 DOF linear computer 


program. 
For pitch motion, а; - coefficients are 
а) = DTHZ 
а, = ОТНТН 
agg = DTHDTH 
ал; = DTHMB 


2. Plenum Air Flow Sensitivity Coefficients 
These coefficients are calculated directly from [2]. 


The ajj S related to plenum air sensitivity coefficients are 
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DMBMB 


f» 
ll 


2 


3. Sway Sensitivity Coefficients 


Before deriving the sensitivity coefficients for 
sway motion, some partial derivatives related to the rudder 
and sidewall hydrodynamics must be calculated. From equation 
(11-11 (Ό)) Y, is already known, this follows the partial 


derivatives of Y, with respect to the state variables 


= r — - 
Ro = эр ο, (64/3) 2 (III-11) 
ЗУ, 
Ry к^. hin 05 C. (C/U) (III-12) 
oY 
R. = = = - C, (C,/u) X, (III-13) 


From equation (II-39), the sidewall cross flow drag 
terms D. partials are derived in the following equations. 


Let us define first some coefficients 


lic 


IID> 


2 ¿ (0) + sign(V,) b. $ (0) 


The partials are 


9р 2 


5, = 55 = - а, sign(V,) MS (0) (III-14) 
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E = - d; V, ^ (0) b (III-15) 
D, - 2d, d, sign(V,) V,(0) (2, - 2,(0)/2) (111-16) 
D, - -24d,d, У (0) (III-17) 
р. = - 24, 4, sign(v,) V,(0) X, (III-18) 


and d. are defined as 


d 3 


2 


г> 


2А b u (III-19) 


Ne 


A u (III-20) 


where u is the craft speed assumed to be constant. 


For sway mode, the state matrix elements are 


= 1020) = Шр, sin 8) (TIT-21) 


Ó 


where D, is given from equation (III-15). 


ф 


(SX) 


397 p 


= нін 


= ES + (в, %4;2,)віп В + 4.2.1 (ІІІ-22) 
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1 (ЗҮ 
m ον 


299 
E: = a a 
= ES 4 (D, d,)sin B dz) (III-23) 
a = 
911 m or 
= Ев + (р - ах )ѕіпв- ах -πι] 
m Y T 3 s 3 s 
(III-24) 


4. Roll and Yaw Sensitivity Coefficients 

Recalling equations (III-6) and (III-7), the net 
roll and yaw moments appear in both equations of motion. In 
this case, to determine sensitivity coefficients is not easy, 
since the yaw and roll partials are involved in the dynamic 
equation with different inertial terms and coefficients. 

In Sys to simplify the problem, the roll and yaw 
partial derivatives are calculated individually, then combining 
these coefficients and dividing by corresponding inertial 
terms give the system sensitivity coefficients for the roll 
and yaw modes. 

а. Roll Mode 

The net roll partial derivatives are defined by 
k.. variables. 


17 
The partial derivative with respect to z is 


πας 29] 


k = — = yy (Ko + K τ T Ка) 


In equation (III-26) all roll components have 


already been derived in the preceding sections. 
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From 


partials are 





equations (II-30) and (II-31), the seal 


9K £ bo 3 | 

—— = - (35) p g(1/sin 31°) Q (0) 

sa = - CN 0 g(1/sin 32?) ф (0) 

For buoyancy partials, let us define the following 
coefficients; 

ο ο οσο (їїї-26) 
b, = 204 Le πι ο) 
p, Š ΝΗ; + 5: (111-28) 
b, à c * b. (III-29) 
be A з“ b. (III-30) 


From equations (II-35) and (II-36), the buoyancy 


partials are 


Ку = 
ΟΖ 


ΟΚ 
97 


DES 


= - z— (bs +b.) 910) 


tan By 


> 81 
= - gan E Pe +249 (0) 
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Equation (II-26) follows; 
aK 2 ӘР, 
552 = -[A (0) P) (0) + (Z. -& (0) =] Ф (0) 


Similarly for sidewall cross flow drag partial is 





οκ: 
s ¿= D, (5. cos B - Z, sin 8) 
The roll partial derivatives with respect to 
8 are 
_ OK 
ЕНТ 155 
y d: = 
= 356(К5ғ + Esa + Kpg + Ка! en) 
where 
ЭК р 
sf _ Wan 3 Š 
πο. = (3734) e g(%,/sin 31%) %(0) 
ӘК b 
Sa ELS w ,3 : 
oK b. b 
DL “2 
mE na зБ 
oK b, b 
ba _ | 2 4 
56 Бап 8-96 7 54 43) 900) 
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15 


where 


75 
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The roll 


The roll 


partial derivative with respect to М, 


ЭК 
--Р 
ο 
ӘР, 
(0) (2. -2.(0)) --т- ф (0) (ІІІ-32) 
5 а ΟΜ. 


partials with respect to $ are 


36 Kor = Esa р Кы Е Кра Ы K. i E 


(111-353) 


3 


(b. /2-3) ° pe g x 


P, (0) (2. = 14 09) g b, 


5S3 





where 


where 


озу = (D, 7 d4) (bo cos fg - 2 sin ß) + d, X. Z 


dK 
_ а _ D 


$ $ 


(б cos ß - 2 sin ß) 


The roll partials with respect to p are 


_ 9Κ 
КИШ “зр 
OE - 
- Den + K, + K,) (III-34) 
o K 
Sw _ _ E 2 5 2 
ας, d, 2 Ay, u 32/1 5. Вс 
oK; зү, 
"o - (Do + =s) (bo cos B - 2 51п 8) 


The roll partials with respect to V are 


A 

II 
co| o> 
als 


79 


(K + K 


en a + К.) (III-35) 


ze 


S 
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524 
r 
ον 


Cy (C/U) 7 


r 


The roll partials with respect to yaw rate are 


711 


where 


ӘК 
SW 


or 


dK 
SW 


or 





ӘК 


or 


b. Yaw Mode 


C. (C. /u) x. Z 


эк 
or 

RR + K, +K) (III-36) 
or sw d r 

а; 2T a 


(5. - 4-Х_) (6. cos B - 72 біп В) 


E 


The following coefficients are defined for yaw 


rotational motion; 


ni £ 2 
2 

2 
па = 


22 4(1 - A) (III-37) 
- T/5 pu & (III-38) 
- sign(V,) о Са 2 ¿ (0) (III-39) 
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Ne 
ω 


σι 2 51 £ 

À 2 
с. = 5. Q 

A 2 
Xx = 4/2 76 Ce u 4 


(III-40) 


(III-41) 


(II1-42) 


The net yaw partial derivatives are defined by 


n.. variables. 





1j 
The yaw partials with respect to q are 
oN 
= SW 
"116 $ 
oN 
x l. 3 z 
96 Ç 5ф pw X ow! b 
where 


IN 

TE 2 
noms 4 n, V (0) b. $ (0) 
None XN) eds ах 
930 ` pw sw w 


(III-43) 


(III-44) 


The Mee and N. partial derivatives can be calculated from 


equations (II-44(a)) and (II-11(d)). 
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The yaw moments partial derivatives with respect 


to V and r are derived similarly, the nij coefficients 


related to these modes are 


_ 8N _ 
Πο ^ 5v (111-45) 
_ ƏN к 


с. Roll Sensitivity Coefficients 
The roll sensitivity coefficients corresponding 
to state matrix elements а. 3 are determined from equation 


(III-7). They are 


ак 

843407 з 

a Kae Ty 

S MENS A nA 
Алл “ОЕ Б nu Е 
ат = Krıık т 


where inertial terms I and I. are found from equations (III- 


8(a)) and (III-8(c)). 
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d. Yaw Sensitivity Coefficients 
For yaw motion, the sensitivity coefficients 
Er 5 have been determined from equation (III-6). 
These coefficients are listed below; 


Sagi, = Si 

ΙΙ, τ 

ars курт 

τε ылу 
τη 
nn cnc oe 


where the inertial term I, is calculated from equation 


ELII-8 (b)). 


D. STATE SPACE REPRESENTATION AND OUTPUT EQUATION 


The state vector of the linear 5 DOF model is shown as 
х = [dz dW d6 dq dM, d$ dp dy ἂν ἂψ dr]? (III-47) 
Then, the system state equation must be 
x = Ах + Bu | ` (III-48) 
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And the output equation is 


he 

ll 
OQ 
>< 


(III-49) 


"A" matrix elements have already been determined. However 
for "B" matrix elements, the step input to the system can 
be considered as a weight removal disturbance and for the 
roll and yaw modes, a rudder order disturbance. 


For weight removal forcing function, the input is to be 


We = W. 
5-1 =‹ Ж a (III-50) 
where 
We = final craft weight 
W; = initial craft weight 
m = craft mass corresponding to the final 
weight We. 


When the rudder order is used, sway, roll and yaw forcing 


functions are generated as follows; 


7771 ! 2 
мау bo) = = (2 X, sin ὃ, + C. 8.) (ІІІ-51) 
roll b71 = (2 X. sin Š, + C. dl = Z / 1, + x /I) 
yaw b111 = (2 X. sin б +С, $5) (x, /I, + 2/1) 
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Finally, defining state variables ху, Kar see, etc., 
the state equations of the 5 DOF linear model are given 


as follows; 


ху = X, 

S ma 22454 | 22575 

ху = х, 

хы 5 аджу % а,;х; % а,,х, % а,5х; 

хы 5 ах) % а;4% 

xc = x, 

ху = 27,X] + 273%3 + 275X5 + 27¿X¿ + 277%7 
πο ο aS 

хз = хо 

Xo = AggXg + Ag7X] + AggXy + Agy1X]1 

X19 ^ Xp 
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ла + + 


а атт οι “Ίος 


"Tod m T 


m ЧОЕ СҮҮС Ө ttle TL 
All coefficients of the state equations are defined 
in the preceding sections. 


Considering outputs of the linear model, only incremental 


changes of the variables can be observed. 
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IV. 


LINEAR FIVE DEGREE-OF-FREEDOM 
COMPUTER PROGRAM DESCRIPTIONS 


The Digital Simulation Language (DSL/360) is used for 


the linear 5 DOF XR-3 computer program. Eleven first order 


integrator blocks are needed for five degree of freedom simu- 


lations and for all linear runs, fourth order Runge-Kutta 


with fixed interval integration method is selected. 


In addition to the craft dimensional values, initial and 


final steady state operating points are entered into the 


program for linearization purposes and plot outputs. 


The input parameters of this program are 


V = 


PHI = 


FYP = 


VL = 


Craft speed (ft/sec) 

Final craft weight (lb) 

Initial craft weight (1b) 

=  Rudder input (positive for port turn) 
Average deadrise angle 


Final steady state value of plenum 
pressure (lb/ft ) 


Final steady state craft draft (inches) 
Final steady state roll angle 


The propulsion system thrust corresponding 
to the craft speed (for each of propellers) 


Steady state lateral speed (ft/sec) 


Steady state yaw rate (rad/sec). 
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The steady state values of the variables are determined 


from 6 DOF nonlinear model simulations. 
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V. SIX DEGREE-OF-FREEDOM NONLINEAR MODEL SIMULATIONS 


In this chapter 6 DOF nonlinear model simulations are 
discussed. 

At first,some of the runs are made in order to determine 
craft initial steady state operating point corresponding to 
craft speed and thrust inputs. Since 6 DOF model weight 
transient response can be analyzed by means of entering proper 
craft initial values into the program, this is simply an 


initial condition response of the nonlinear dynamic system. 


A. WEIGHT REMOVAL TRANSIENT 
For 30 and 20 knot craft speeds, the steady state values 


of some variables are listed in Table I. 


Table I 
Steady State Values of the Weight Removal Response 


P 


Speed Pitch b Draft Thrust Weight 

(Кое) (degree) (1Ь/ҒЕ2) (inch) (1b) (1b) 
30 0.95 24.8 8.6 355.6 6722 
30 0.26 24.8 5.6 287.8 6050 
20 0.45 24.8 6.06 186.2 6050 


For 30 knot craft speed, two different steady state pitch 
angles are observed corresponding to the different weights. 
In [2] for pitch transient response, 0.65° steady state pitch 


angle value has been used for both 11пеаг and nonlinear 
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simulations. On the other hand as seen in Table I, 0.26° 

and 0.95° steady state pitch values are obtained. This dis- 
agreement may be explained by noting that different dimensional 
values were used or more probably a different modified 6 

DOF computer program deck was used. 

For 6050 lb craft weight, the steady state pitch angle 
0.26° is in agreement with the value determined by Menzel in 
figure 6 in [6]. 

During the weight removal response, simulations of 6 DOF, 
two different stern seal base leakage areas are used. In [1] 

a smaller leakage area is assumed for the linear heave only 
model as А, = 0.384. However А, = 3.84 also has been simulated 
and it is understood that the smaller leakage area is more 
reasonable considering other steady state variables in the 

air flow dynamic equations (II-23) and (II-24). Since at 


the steady state condition Qin = Q must be satisfied, that 


out 


permits only smaller A, for a given Q, and p steady state 


values. Therefore А, = 0.384 is preferred and used for 


nonlinear simulations. 


B. TURNING MANEUVER 

When the 3 DOF pitch, heave and roll linear model was 
developed, the 6 DOF nonlinear model simulations were used 
in order to determine all steady state roll moments during 
the turn maneuver. In addition to the. roll moments, yaw 
and sway outputs were also investigated and the results are 


tabulated in Tables II, III and IV. 
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Table II 


Steady State Roll Moments 


(rudder = 


Roll Moments 
Stern seal FKSS 
Bow seal FKBS 
Sidewall ЕКЅИ 
Rudder  FKRUD 
Propulsion  FKP 
Aerodynamic  FKAED 


Plenum  ABPB*PHI*(-Z2) 


59, 30 knot speed) 


Positive 


131; 


0.057 


215.3 


Table III 


Steady State Yaw Moments 


(rudder 


Yaw Moments 
Stern seal FNBS 
Bow seal  FNSS 
Sidewall FNSW 
Rudder FNRUD 
Propulsion FNP 


Aerodynamic  FNAED 


= 5°, 30 knot speed) 


Positive 


3.22 


0. 


532.65 
2-18 
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Negative 


=22,1 


Negative 


0. 
-482.2 


- 56.41 





Table IV 
Steady State Sway Forces 


(rudder = 5°, 30 knot speed) 


Sway Forces Positive Negative 
Sidewall FYSW _ =351.8 
Rudder  FYRUD - - 47.6 
Aerodynamic FYAED - - 19.07 
Centrifugal  -R*U*AM 419.24 = 


From Table II, it is understood that for a given (+) 
rudder input, the sign convention of the craft dynamics was 
in agreement with turning to the left. However, a negative 
sway displacement plot is observed for this maneuver. As a 
result, in the R.H.S. subroutine of the 6 DOF program, the 
Sign correction should be done for y-displacement. 

The steady state sway forces acting on the XR-3 test 
craft in a turn to the left are shown in figure 6. As seen 
in this figure and from Table IV, the rudder sway force is 
in the steady state condition. For a left turn, initially 
the rudder sway force is generated in the reverse direction 
and finally a negative steady state rudder sway force is 
produced. Further investigation is done in the RUDDER sub- 
routine dynamic equations of the 6 DOF computer model, but 
no error has been determined in these equations, Also, 
different craft speeds and different effective rudder angle 
coefficients (ENDFAC) are used but everytime the same 


configuration is observed. 
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Steady state sway forces configuration 
at the left turn 
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In addition to rudder dynamics, the propulsion system 
of the XR-3 contributes to the roll, yaw and sway motions. 
However this contribution is not a dominant factor compared 
with the rudder effect. 

For the XR-3 test craft propulsion system configuration, 
the direction of the thrust vector is also changed as well 
as the rudder deflection. As seen from Tables II, III and 
the simulations with different rudder order, it is understood 
that the propulsion thrust component due to the rudder order 
is not calculated properly, therefore the following modifica- 


tion is made in the PROPULSION subroutine of the 6 DOF program, 


Күс DSSTUISTSSCBILTHSS*GCD 


where FYS represents the sway force due to the starboard 
propeller, and the others are defined in the propulsion 


subroutine program. 


C. DEADRISE ANGLE EFFECT 
As reported in [4], the deadrise angle has a major effect 
on the surface effect ship's roll behavior and therefore is 
an important design consideration. The deadrise angle effect 
on the roll motion has been investigated by Reidel in [3]. 
In this investigation, some arbitrary coefficients (PROMO 1, 
2, etc.) are used for the vertical force component of the 
sidewall hydrodynamics. From this investigation, the conclusion 
reported by Reidel was that more damping could be obtained for 


the pitch and roll motions by changing the deadrise angle. 
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The deadrise effect is reinvestigated by simulating the 
system with different rudder orders including zero disturbance, 
then the heave forces due to sidewall hydrodynamics were 
compared to each others. As a result, no significant sidewall 
heave force change was observed. For instance, for 0” and 
5° rudder inputs, the sidewall heave forces (FZSW) are found 
to be -568.8 and -568.05 respectively. This result is not 
acceptable, since in a turn maneuver, the cross flow drag 
is generated at the outboard sidewall, so the deadrise projec- 
tion of cross flow drag must be contributed to the sidewall 
heave force. More investigation has been made in the SIDEWALL 
subroutine of the 6 DOF program and it has been determined 
that only the force component due to added mass given by 
equation (II-21) is projected as a lift force. Actually the 
cross flow drag component must be included to this projection, 
since the frontal area for cross flow drag component is out- 
board of the submerged section of the sidewall. Therefore 
it is concluded that the XR-3 test craft 6 DOF program SIDEWALL 
subroutine does not include effectively the lift force due 
to deadrise. On the other hand, the lift force and horizontal 
drag component can have different drag coefficients, therefore 
the modification of the sidewall subroutine is left for future 


Studies. 
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VI. COMPARISON OF NONLINEAR AND LINEAR MODEL RESULTS 


A. C.G. ACCELERATION (see figs. 9 and 23) 

A 10% weight removal disturbance is used for both models. 
The linear and nonlinear model showed very rapid response 
with small rise time. As seen in figure 9, the 6 DOF model 
transient is relatively more damped compared with the linear 
model response. The linear model showed a small overshoot 
but it is also sufficiently damped. Small oscillatory behavior 
is observed in both figures which is due to pitch-heave 
coupling. In general the two models agree in transient response 
and steady state behavior. 

Both of the C.G. acceleration curves were normalized by a 


factor "g",the gravitational constant. 


B. PITCH TRANSIENT (see figs. 10 and 24) 

As seen in these figures, both models showed an oscillatory 
transient response with the same frequency, w = 3.12 rad/sec. 
However their steady state values are different. In fact, 
the pitch moment due to plenum pressure and the effective 
plenum pressure center play a dominant role for pitch motion. 
In [2] Barnes developed plenum modeling including effective 
plenum roof area calculations and plenum pressure centroid 
shift due to draft variation. Therefore the linear model 
plenum modeling is different from the nonlinear modeling, 
so that a greater plenum moment arm is calculated in the linear 


model. 
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The steady state pitch angle values are 0.5° and 0.75° 


for nonlinear and linear models correspondingly. 


C. PLENUM PRESSURE TRANSIENT (see figs. 7 and 21) 

As seen in these figures, both models are exactly in 
agreement with each other. In each curve, the plenum pressure 
suddenly drops to 21.6 (1b/£t^) then increases up to its 
Steady state value with sufficient damping. This is expected 
behavior because the plenum volume increases when 10% weight 
is removed, obviously this disturbance causes a sudden drop in 
plenum pressure. As seen before, pitch-plenum pressure coupling 
had been observed during the exponential increasing phase of 
the plenum pressure. At the fifth second of simulation time, 
nonlinear model plenum pressure is about 24.29 while the 


linear model pressure is 24.25. 


D. DRAFT TRANSIENT (see figs. 8 and 22) ` 

Excellent agreement was achieved between linear and non- 
linear model draft responses. As expected, the draft decreases 
exponentially when the weight is removed from the craft. 
The draft is found to be 5.7 inches for the linear model, 
5.85 inches for the nonlinear model at the fifth second of 
simulation time. 

As a result, both models exhibited the same draft 


responses. 


E. ROLL RESPONSE 
As mentioned before, at first a 3 DOF linear model (pitch- 


heave-roll) had been developed. 3 DOF linear and 6 DOF 
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nonlinear models roll behavior comparison is done by means 
of a step input roll moment disturbance to the linear model 
and a turning command to the nonlinear model. In figure 13 
6 DOF nonlinear model roll angle response is shown for a 5? 
rudder order (left turn),in this simulation 1.18? outboard 
Steady state roll angle is generated. From the R.H.S. sub- 
routine of the 6 DOF program, all the steady state roll moment 
components are determined as seen in Table II and the sum of 
these disturbance moments applied to the linear model roll 
mode as a step input. Roll behavior is analyzed from figure 
27. As expected, the linear model showed a greater overshoot 
and more oscillatory response compared with the nonlinear 
roll transient. Since the linear model has been disturbed 
by a step roll moment corresponding to the nonlinear model 
steady state roll moments. On the other hand, its steady 
State roll angle is found to be 1.28?. The roll steady state 
difference of two models is about 8.4%. For a 3? rudder input, 
nonlinear and linear simulations are obtained in figures ll 
and 25 with the steady state values 0.62? and 0.65? respec- 
tively. Hence the difference is 4.8% in the steady state 
values for these simulations. 

Figures 12 and 14 show XR-3 yaw rate plots and corres- 
ponding linear simulations are shown in figures 26 and 28. 

As a result, the 3 DOF linear model and 6 DOF program 
roll simulations are in agreement with each other within 


10% difference. 
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As discussed in the previous chapter, some problems 
are observed in the 6 DOF nonlinear program. Therefore 
developing two more modes for the 3 DOF model, the 5 DOF 
linear model has been generated. By means of this model, 
roll behavior of the XR-3 in a turn, deadrise angle and 
cross product inertial term effects have been investigated 
for 20 knot craft speed with different rudder deflections. 

The roll response can be seen in figure 29 omitting the 
cross product inertial term I. = -2800.) in the equations 
of motion, applying a 10° rudder order and assuming 70 degree 
average deadrise angle of the craft. For this particular 
test, the craft starts inboard rolling with the oscillatory 
transient then goes to 0.47° steady state roll. Using the 
same input parameters and also including the cross product 


term I, in the equations of motion, smaller overshoot (18%) 


2 
and a more damped roll response has been observed in figure 
31. As seen in this figure, at the initial phase, secondary 
overshoot has been generated. The steady state roll angle is 
0.40° and the oscillation frequency w = 3.5 rad/sec., for 
this simulation. On the other hand, the corresponding 6 DOF 
nonlinear simulation is shown in figure 17 and also a rc11 
rate plot is shown in figure 18. The nonlinear simulation 

in figure 17 shows a more damped transient but greater steady 


state roll behavior (1.32°) compared with the corresponding 


linear simulation shown in figure 31. 
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Figures 30 and 32 present linear simulation curves for 
5 and 15 degree rudder inputs respectively. Both curves 
exhibit the same linear model response discussed before. 
Finally the linear and nonlinear model steady state responses 


are summarized as follows. 


6 DOF 5 DOF 
Rudder Nonlinear Linear Reference 
5 0.57 0.24 figures 15 and 30 
10 1332 0.40 figures 17 and 31 
15 2.2 0.45 figures 19 and 32 


The deadrise angle effect is investigated by using differ- 
ent deadrise angle values for 10 degree rudder input and 20 


knot craft speed simulations. The results are tabulated in 


Table V. 
Table V 
Sidewall Deadrise Angle Effect 
(rudder = 10°, 20 knot speed) 
Lateral 
Deadrise Roll angle Speed Yaw rate 
angle (degree) (degree) (ft/sec) (rad/sec) 
45 -1.01 4.18 -0.064 
50 -0.84 4.04 -0.069 
59 -0.62 3.90 -0.076 
60 -0.30 2283 -0.080 
65 0.02 3.73 -0.085 


70 0.40 3.60 -0.100 
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For the linear model, the lateral velocity and yaw rate 
plots are shown in figures 33 and 34. Both of them showed 
sufficiently damped transient response. From figure 33, 
the lateral speed steady state value is found 3.6 (ft/sec) 
and as expected for turning to the left, the negative yaw 


rate is generated as -0.1 (rad/sec). 
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Figure 7. 6 DOF Plenum pressure transient 
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Figure 8. 6 DOF draft transient 
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Figure 9. 6 DOF C.G. Acceleration 
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Figure 10. 6 DOF pitch angle transient 
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Figure 12. 6 DOF roll rate plot, oe = 3° 
30 knot speed 
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Figure 13. 6 DOF roll angle plot, ó, = 5° 
30 knot speed 
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Figure 14. 6 DOF roll rate plot, δ, = 5° 
30 knot speed 
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Figure 15. 6 DOF roll angle plot, δ. = 5° 
20 knot speed us 
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Figure 16. 6 DOF roll rate plot, $, = 5° 
20 knot speed 
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Figure 17. 6 DOF roll angle plot, ó, = 10° 
20 knot speed 
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Figure 19. 6 DOF roll angle plot, S = 15° 
20 knot speed 
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Figure 20. 6 DOF roll rate plot, ô = 15° 
20 knot speed 
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Figure 21. Linear 5 DOF plenum pressure transient 
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Figure 22. Linear 5 DOF draft transient 
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Figure 24. Linear 5 DOF pitch angle transient 
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Figure 25. Linear 3 DOF roll angle plot, ο. = 3° 
30 knot speed 
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Figure 26. Linear 3 DOF roll rate plot, δ. = 3° 
30 knot speed 
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Figure 27. Linear 3 DOF roll angle Do, 6 = 5° 
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Figure 30. Linear 5 DOF roll angle plot, 6, = 5° 
20 knot speed 
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Figure 3l. Linear 5 DOF roll angle doc δ. = 10° 
20 knot speed 
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Figure 32. Linear 5 DOF roll angle plot, 6. = 15° 
20 knot speed Е 
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Figure 33. Linear 5 DOF lateral velocity plot 
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Figure 34. Linear 5 DOF yaw rate plot 
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VII. CONCLUSIONS AND RECOMMENDATIONS 


In the preceding sections, the 5 DOF linear XR-3 model 
was developed, then weight removal response and roll behavior 
are compared with 6 DOF nonlinear simulations. Also some 
of the craft sensitivity parameters are investigated for the 
XR-3 test craft. 

The results of this study are summarized and recommenda- 
tions are made as follows. 

1. The system response is investigated by both a linear 
model and a 6 DOF nonlinear model by means of sudden weight 
removal from the C.G. and by sway force disturbance in a 
turn. 

Comparison of results shows that the difference is 
within 10% for both models. In fact, assuming accuracy in 
measurements of about 10%, the linearized model gives 
acceptable results as compared with the nonlinear model. 
Obviously the linear model simulations take less computer 
CPU time. 

2. Concerning roll response during the turn maneuver, 
the product inertia term (Izz) affects transient response 
such as overshoot, damping, etc. More overshoot and 
oscillatory response have been observed without Te: 
However, this term is found to be -2800. from the 6 DOF 


program. It is suggested that this number be verified. 
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3. The importance of deadrise is observed for roll 
steady state behavior. The XR-3 test craft steady state 
roll angle depends upon the average deadrise angle of the 
sidewall. Different average deadrise angles have been 
simulated by means of the linear model and it is observed 
that the craft rolls inboard or outboard according to the 
sidewall deadrise configuration. 

The 6 DOF nonlinear model SIDEWALL SUBROUTINE 
must be modified by including a cross flow drag component 
deadrise projection (lift force). However a new drag 
coefficient is needed for this lift component. 

4. The aerodynamic effect of the craft contributes 
about 5% of the total yaw moment. 

5. The XR-3 failure mode operation (such as fan failure) 
can be simulated by using either the linear or nonlinear 
model. But, the nonlinear model is preferable for failure 
mode simulations, since it includes sidewall and stern 


seal leakage configurations. 
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